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Schematic  diagram  of  an  artificial  Kerr  medium. 
Particles  at  the  focus  are  moved  by  the  optical 
dipole  force  to  change  the  refractive  index. 

Schematic  diagram  showing  the  optics  for 
calculating  the  dynamic  nonlinear  coefficient  of 
an  AKM.  (From  Ref.  7.) 
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The  estimated  in  vs  t  curves 
laser  focused  at  AO- urn  diameter  a)  In  a  llqi 
mixture  and  b)  in  an  AKM.  c)  A  corresponding 
curve  for  a  I  W/cm^  power  density  In  the  liquid 
mixture.  Specifications  of  the  liquid  mixture 
and  the  AKM  are  given  in  the  text. 

Optical  configuration  for  calculating  the  focal 
length  of  a  thermal  lens  and  the  transmission 
through  aperture  A. 

Dynamic  nonlinear  transmission  through  aperture 
A  in  Fig.  A  due  to  thermally  Induced  self- 
defocusing  for  A  -  0.055  (upper  curve)  and  0.002 
(lower  curve)  cm.  Absorption  loss  in  the  liquid 
is  1AZ;  overall  transmittances  in  the  absence 
of  a  laser  (l.e.,  t  -  0)  are  85  and  7AZ, 
respectively. 

Schematic  diagram  of  the  nonlinear  interface. 
Media  a  and  b  occupy  the  regions  Z  >  0  abd  <  0, 
respectively. 
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Note  that 


for  various  incident  angles  9^. 
decreases  as  I  increases. 


FOREWORD 


This  document  is  the  first  volume  of  the  Final  Report  (two  volumes) 
of  Contract  No.  DAAK70-84-C-0022,  "IR  Artificial  Kerr  Media,"  from  U.S. 
Electronics  R&D  Command,  Night  Vision  and  Electro-Optics  Laboratory 
(NV&EOL).  The  work  was  performed  by  W.P.  Chen  and  C.C.  Frazier  with  the 
assistance  of  H.  Jones,  M.A.  Harvey,  L.M.  Gillette,  M. P.  Cockerham,  and 
P.L.  Porter  during  the  period  18  April  1984  to  31  July  1986  under  the 
management  of  J.M.  Chen. 

The  authors  are  grateful  to  G.  Wood,  E.  Sharp,  and  R.  Shurtz*  of 
NV&EOL  for  their  support,  encouragement,  and  fruitful  discussion.  They 
also  wish  to  acknowledge  M.  Kroll  for  pointing  out  the  thermal  lenslng 
phenomenon  and  S.  Guha  for  elucidating  it. 


I.  INTRODUCTION 


In  recent  years ,  nonlinear  optical  materials  having  large  nonlinear 
optical  coefficients  when  combined  with  novel  optical  mechanisms  have  found 
many  applications  in  the  field  of  phase  conjugation,  image  processing, 

®  optical  computing,  and  optical  switching. In  this  volume  of  the  Final 

Report  we  will  discuss  two  types  of  nonlinear  materials,  artificial  Kerr 
medium  (AKM)  and  liquid  mixture;  and  two  types  of  nonlinear  optical 
mechanisms,  self-defocusing  and  nonlinear  interface. 


e 


i 
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II.  NONLINEAR  MATERIALS 

A.  ARTIFICIAL  KERR  MEDIUM 

An  artificial  Kerr  medium  is  a  suspension  of  transparent  particles 
in  a  transparent  liquid  of  a  different  refractive  index. a  focused 
laser  beam  can  exert  substantial  radiation  pressure  on  the  particles  and 
cause  nonlinear  effects.  Under  the  influence  of  such  pressure,  the 
particles  redistribute  (Fig.  1). 

When  a  dielectric  particle  is  in  an  electrical  field  £,  it  behaves 
like  a  dipole  with  a  polarization  P  given  by^^ 

?  -  p£  (1) 


where  p  is  the  effective  polarizability  of  the  particle.  For  a  sphere,  the 
p  is  given  by 


2  \  3 
P  -  (__)  a 


(2) 


r 

where  nr  *  na/n£,  *s  tbe  ratio  of  the  refractive  indices  of  the  sphere 
and  the  liquid  (n^)  and  a  is  the  radius  of  the  sphere.  The  sphere 
experiences  a  force  due  to  the  optical  field  gradient  given  by^^ 


(ng) 


and 


£  -  (?  .  7)  £  -  (1/2)  p  t£q2  =  -V<p 


i 


(3) 

(4) 


where  EQ  is  the  time-averaged  rms  electric  field  of  the  light  wave  and  4>  is 
the  light-induced  potential. 

For  a  light  beam  incident  onto  a  liquid  suspension  of  spheres,  the 
optically  induced  force  moves  the  spheres  with  ng  >  n£  and  ng  <  n^  into  and 
out  of  the  high-optical-intensity  region,  respectively.  In  both  cases,  the 


motion  of  the  spheres  creates  refractive  index  changes  which  are  analogous 
to  those  in  a  medium  with  a  positive  optical  Kerr  coefficient,  and  the 
effective  refractive  index  for  the  suspension  system  is  given  to  a  good 
approximation,  as  long  as  f  «  1,  by 


n  ■  n .  +  f  (n 

medium  £  s 


(5) 


where  f  (»  — a^N)  is  the  volume  fraction  of  the  spheres,  and  N  is  the 
number  of  spheres  per  unit  volume.  For  f  <  0.1,  the  given  by  Eq. 

(5)  will  deviate  from  that  derived  from  the  more  accurate  Maxwell-Garnet 
theory  by  no  more  than  2%. 


For  application  of  this  medium  as  an  optical  switch,  the  scattering 
loss  associated  with  the  nonlinear  effect  must  also  be  considered.  When 
the  radius  of  the  spheres  is  much  less  than  the  wavelength  of  the  incident 
light,  the  scattering  loss  per  unit  length,  uQ,  is  given  by  the  Rayleigh 
formula^) 


(6) 


where  XQ  is  the  wavelength  of  the  incident  light  in  vacuum.  One  figure  of 
merit  of  the  medium  is  the  ratio  of  its  nonlinear  coefficient  to  Its 
scattering  loss.  This  will  be  discussed  in  the  following  sections  for  the 
static  and  dynamic  cases.  We  can  rewrite  Eq.  6  as 


a  X 
o  o 


(7) 


As  we  keep  f  constant  and  scale  a  with  XQ,  the  aQX0  can  be  kept  constant. 


1 .  Continuous-Wave  (Static)  Case 


In  the  static  case,  the  flow  of  the  spheres  induced  by  the  optical 
field  is  counterbalanced  by  thermal  diffusion.  At  equilibrium,  the  ratio 
of  N  inside  and  outside  the  beam  is  proportional  to  exp  (-$/kT),  where  $  is 
the  optically  induced  potential  inside  the  beam  relative  to  that  outside 


* 


and  is  proportional  to  the  optical  power  density  I  of  the  light  beam. 
Usually,  p  «  kT,  and  we  can  expand  exp  (-^/kT)  as  ( 1— 4>/kT) .  Thus,  the 
optically  induced  difference  in  the  sphere  density  is  ~  N4>/kT ,  and 


n  aa  m  no  +  4~  (n  ~  n  )  -  4“  (n  -  n  )  . 

medium  £  3  s  13  s  £  kT 


(8) 


Because  4>  oc  I,  the  third  term  in  the  equation  describes  the  effective 
optical  Kerr  effect.  The  static  effective  optical  Kerr  coefficient  rt2  is 


then  given  by 


(6) 


n . 


n 


2  "nr 


(nr2  -  1)  a3  f 


2  ckT  (nr  +  1)  (n2r  +  2) 


(9) 


When  the  n2  is  normalized  to  the  aQ  (scattering  loss)  and  the  aQ \Q 


to  derive  the  figure  of  merit,  we  obtain: 


o 

n. 


16ir3ckTn  2 


n  +  2 
(— r — : — r) 


a  X 
o  o 


n  +1 

*  3  r  n  2  +  2 
°  — . — > 


16ir3ckTn  2 
£ 


n  +  1 
r 


( 10a) 

(  xOb) 


Note  that  the  proportionalities  ^2^  ao  an(^  n2'/(Xo^o  vary  38  the  fourth  and 
third  power  of  the  wavelength,  respectively.  As  long  as  the  scattering 
loss  is  kept  constant,  n2  Is  independent  of  a,  and  the  redistribution  of 
spheres  is  determined  by  the  thermal  equilibrium.  Therefore,  the  resultant 
static  n2  depends  on  temperature  and  the  optically  induced  potential 
difference,  and  not  on  the  field  pattern  of  a  focused  light  beam. 


2. 


Pulse  Case 


The  field  pattern  does  play  an  important  role  in  the  pulse  case, 
because  the  steepness  of  the  field  gradient  determines  the  forces  on  the 


spheres  and  their  accelerations.  When  a  TEMQo  laser  beam  is  focused  as 


shown  in  Fig.  2,  the  field  pattern  at  any  plane  z,  expressed  in  Gaussian 
units,  is 


-5- 


i 


/8tt  p *  \  1  / 2 
vcn  7 


ux 


Z  (r,z)  -  exP  H  Ik*  "  ♦(*)]}  exp 

At  the  focal  spot  (z-0),  v  (r,o)  is 


(11) 

(12) 


UQ  (D 


/li- 


,  .  k  r% 
x  w  exp  ("J  2  " > 
0  qo 


(13) 


where  P'  ■  total  input  power 

n  -  refractive  index  of  the  medium  at  low-power  irradiation 
c  -  speed  of  light 

wQ  -  radius  of  the  beam  at  the  focal  spot 
k  -  2ir/ X 

X  -  Vn 

q(z)  -  q  +  z 

q  «  1  xw  /X  =  1z_ 

J  o  J  R  _ 

zR  -  the  Rayleigh  length;  at  z  -  ±  zR,  w  »  /  2  wq 

_ 

q  (z)  xw2(z) 


R(z) 


z  + 


[i  ♦  £> 


w(z)  •  u  /  1  +  (|— ) 

°  ZR 

i|»(z)  -  tan”1  (•§-) 

ZR 

//  |  "uo  | 2  2*rdr  -  1. 

The  field  pattern  in  the  region  z  <  0  is  a  mirror  image  of  that  for 
z  >  0  at  the  plane  z*0.  When  a  laser  beam  with  beam  diameter  D  is  focused 
by  a  lens  with  focal  length  L,  the  wQ  is  given  by 


-7- 


c 


(14) 


w  -  — — - 
o  itD 


provided  — —  <<  1. 

TTU 


Now,  from  Eq.  (3),  we  can  calculate  force  as 


?  -  ~  p  VE  2 
2  K  o 


-2pr  |E 


.  1  - 


[‘ *  t)2J 


mil  >Eo’2 


Along  the  optical  axis  (r-0),  the  force  Is  only  In  the  z  direction  and 
towards  the  focal  point,  with  maxima  occurring  at  |z|  »  — zR.  At  the 
focal  plane  (z-0),  the  force  is  only  on  the  plane  and  towards  the  focal 
point,  and  a  maximum  occurs  at  r  -  0.5  wQ.  In  a  situation  where  a  sample 
Is  thinner  than  zR  and  locates  at  z-0,  the  force  is  mainly  in  the  radial 
direction. 


Substituting  p  from  Eq.  (2)  into  the  first  term  of  Eq.  (15),  we  can 
explicitly  derive  the  force  at  z-0  as 


| F |  *  r“»f2 

11  cn  l 


2  ' 
n  -  1 


n  + 
r  „ 


V- Da3J  "Texp  r^T 


=  Aa3_E_exp  p. 

with  the  force  toward  the  center  and  A  equal  to  the  value  in  the  bracket. 
In  cgs  units,  the  equation  of  motion  of  the  spheres  is  given  by 

mr  ■  -6  wa  rtf  -  A  a}  exp  (17) 

where  m  is  the  mass  of  the  spheres  and  n  is  the  coefficient  of  viscosity. 
The  first  terra  on  the  righthand  side  of  Eq.  (17)  is  the  viscosity  force. 


-8- 


In  the  region  r  <  0.5  wQ,  where  the  force  is  significant,  Eq.  (17) 
can  be  approximated  as* 


r 


6iran  ^ 
m 


5r  +  gr. 


(18) 


The  force  due  to  light  is  overestimated  by  Eq.  (18),  so  that  the  actual 
response  will  be  slower  than  the  calculated  one.  In  the  cases  that  we  are 
likely  to  encounter,  namely,  a  <  2  pm  and  wQ  >  3  p,  then  4  6  «  6  ,  and 
the  solution  to  Eq.  (18)  becomes  simply 


r 


-( 8/ S)t 
i 


ri  "-I1' 


(19) 


where  r^  is  the  initial  position  of  the  sphere  at  the  onset  of  the  pulse. 

Now,  we  can  estimate  the  change  of  density  of  the  spheres  in  the 
region  r  <  0.5  wQ  as  a  function  of  time  as  follows.  All  of  the  spheres 
originally  located  within  a  distance  rt  from  the  origin  move  to  within  a 
distance  r^,  given  by 

rf  -  r±  (1-  |  t).  (20) 


Therefore,  the  density  within  r^  increases  by 
(irr  2  -  wr  2)  N 

AN  -  - - - 2 -  -  N  (2  t). 

"f 

The  resulting  change  in  the  refractive  index  is 


(21) 


*  First  pointed  out  by  G.  Wood  at  NVEOL. 


(22) 


V' 


1 


v 


k>i 


•a 


An  =*  AN  (n  -  n  )  -^  ira^ 
S  £  J 


8 it  3  ,  .  8 

—  3  <n8‘  "«>  1  Nt- 


Substituting  5  and  S  from  Eq.  (18)  and  N  from  Eq.  (6)  with  f  -  a^N 


into  Eq.  (22),  we  have 


/  2 
n  z+  2 
r 


6ttV2  \nr+1y  CT1 


-j  t  IaQ  (in  cgs  units) 


where 


1  s-h p’ 

irw 

o 


is  the  average  power  density.  Using  the  relationship  An  *  n^I  and 


rearranging  the  equation,  we  obtain 


3  2 

16u  c  n 

l 


2 

nr  +  2 


nr  +  1 


where  n^  is  the  dynamic  nonlinear  coefficient. 


Like  the  figure  of  merit  of  the  CW  nonlinear  effect,  n2/ao\)  (Eq. 
10),  n'2/ao\>  3  function  of  XQ,  n^,  and  nr>  However,  it  depends  on  the 

viscosity  and  the  beam  width  Instead  of  the  temperature.  It  should  be 


noted  that  the  response  becomes  faster  as  the  sphere  becomes  smaller. 


We  can  illustrate  this  effect  by  calculating  the  responses  for  AKM 


-1-1 


with  the  following  parameters:  n£*  1.5,  ng  -  2.35,  u*  0.02  g  sec-1 cm 

«« J  .  .  ic  1  a  ski. i  g 1  j  _ _ i  .. 


and  aQ  ■  15  cm  .  A  100-pm-thick  layer  of  this  medium  would  cause  only  a 
147  scattering  loss.  With  the  parameters  given,  the  product  a^f  is  then 
fixed  by  Eq.  (6)  for  a  given  wavelength  XQ.  A  small  radius  will  be 
selected  to  produce  fast  response.  However,  since  the  volume  fraction  f 
cannot  exceed  100%,  there  exists  a  lower  bound  of  the  radius  a  which,  for 
this  case,  is  ~  2200  A.  So,  if  we  take  a  ■  2200  A,  then  the  An  vs/t  curve 
given  by  Eq.  (24)  and  shown  in  Fig.  3  sets  the  optimum  response  of  the  AKM 
for  a  0.38-W  CO2  laser  focused  to  a  40-im  diameter. 


B.  LIQUID  MIXTURE 

A  liquid  mixture  can  have  a  very  large  dynamic  nonlinear  effect. 
However,  the  effect  was  not  appreciated  in  connection  with  optical 
switching  limitation  applications.  Absorption  of  incident  light  upon  a 
liquid  can  heat  up  the  liquid  and  change  its  refractive  index.  In  regions 
of  the  optical  spectrum  far  removed  from  the  absorption  bands  of  the  liquid 
and  at  temperatures  away  from  phase  changes  or  critical  points  of  the 
liquid,  the  temperature  coefficient  of  the  index  of  refraction,  dn/dT(K-1), 


dn/dT  -  (3n/3T)  +  (3n/3p)  (3p/3T) 

P 


is  determined  primarily  by  changes  in  the  sample  density  (second  term). 
Most  liquids  expand  when  heated,  resulting  in  a  negative  value  of  dn/dT. 

The  thermally-induced  change  in  n  can  be  estimated  from  the  heat 
produced  by  the  absorption.  The  power,  Ptrans»  transmitted  through  a 
thickness,  £,  (cm)  of  a  sample  with  absorptivity,  a,  (cm-1)  is  given  as 


trans 


where  P  (W)  is  the  incident  power.  For  a£  <<  1,  the  absorbed  power  is 
approximately 


p  .  -  P  -  P„  -  Pa£  . 

abs  trans 


Consider  the  situation  of  a  laser  focusing  to  a  waist  w  upon  a  sample  with 
a  ■  15  cm-1  and  £  -  100  um.  The  absorption  is  ~14%.  The  rate  of  energy 
delivery  to  the  sample,  Qjn,  is  Pa£/J,  where  J  *  4.18  J  cal-1  is  Joule's 
constant.  We  assume  complete  dissipation  of  the  incident  energy  as  heat. 
The  rate  of  temperature  increase,  T,  is 


Pat/Jw2  tpC  -  (j^)  (-iy)  I  j2_I 

p  xw  p 


»  1  *  i  ;  v  ■ .  i  • 


T 


2  3 

depending  on  the  sample  volume  nw  £(cm  ),  the  specific  heat  Cp 
(cal  g-*  K-*),  and  the  sample  density  p(g  cm  3).  The  heat  is  then 
conducted  radially  from  the  illuminated  region  into  its  surroundings.  The 
thermal  conductivity  <  (cal  sec-*  cm-*  K-*)  relates  Qout>  the  time  rate  of 
heat  transfer  by  conduction  to  the  surface  area,  in  the  present  case  2irw)t  , 
to  the  distance  (~  w)  between  the  heated  spot  and  its  surroundings,  and  to 
the  difference  in  temperature  AT;  hence  QQut  m  (<)  (2irwi.)  (AT)/w. 

When  the  sample  temperature  reaches  an  equilibrium  value,  0in  - 
0out*  ^  the  rate  of  temperature  increase  T  calculated  above,  equilibrium 
would  be  achieved  in  a  time 

At  -  AT/T  -  [(  Paf,/ J)w/  le.irwfc]  /(  Pal/  Jttw2  JlpC  )  -V?w2pC  /  <. 

The  expression  for  At  above  represents  a  fundamental  time  for  determining 
whether  the  heat  conduction  shall  or  shall  not  be  taken  into  considera¬ 
tion.  For  a  pulse  shorter  than  At,  we  can  ignore  the  conduction  loss.  Of 
course,  these  calculations  cannot  be  taken  too  literally  and  are  only  meant 
to  provide  an  appreciation  for  the  order  of  magnitude  of  the  effects  being 
considered  here.  Detailed  analysis  can  be  found  in  Ref  8. 

Using  liquid  benzene  with  some  dye  as  a  typical  organic  sample, 

where  l  -  0.88  g  cm--*,  Cp  -  0.41  cal  g-*K-*,  and  ic  ■■  3.41  x  10-^ 

cal  sec-*cm-*K-* , and  assuming  that  dn/dT  in  the  IR  region  is  the  same 

as  that  in  the  visible,  dn/dT  *  -3.9  x  10”^  £-1^(10)  we  can  estimate  the  T 

2 

and  At  for  a  1-W  laser  focused  to  an  area  of  1  cm  .  The  resultant  T  is  ~ 
9.95  K  sec-*.  For  the  case  of  short  laser  pulses  in  which  we  can  ignore 
conduction,  the  T  will  cause  a  decrease  of  n  given  by 

ft  -  -3.9  x  10’3  (W/ cm2 )~* sec”1 .  (29) 

The  change  in  refractive  index,  An,  of  the  sample  is  isotropic.  The  An  vs 
t  curve  (curve  a)  of  the  typical  example  is  also  shown  in  Fig.  3  and  can  be 
compared  with  the  corresponding  AKM  curve.  We  see  that  the  photothermal 
effect  is  at  least  three  orders  of  magnitude  higher  than  the  artificial 


# 


c 
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Figure  3.  The  estimated  An  vs  t  curves  for  a  0.38-W  Cr>2  laser  focused  at 
40-un  diameter  a)  in  a  liquid  mixture  and  b)  in  an  AKM.  c)  A 
corresponding  curve  for  a  1  W/cm^  power  density  in  the  liquid 
mixture.  Specifications  of  the  liquid  mixture  and  the  AKM  are 
given  in  the  text. 


Kerr  effect  at  the  same  response  speed.  Also  shown  is  the  An  vst  curve 
(curve  c)  at  lW/cm^  power  level.  It  is  worth  noting  that  the  photothermal 
effect  (for  t  >  1  msec)  is  comparable  to  the  strong  nonlinear  optical 
effect  in  a  semiconductor  at  the  band  edges.  For  example,  n2  of 
Hg0  gCdp^Te  (having  Ag  *  7.5  um  and  a  ■  26  cm-*  at  10.6  um)  is  8  x  10-® 
(W/cm^)-*  for  180  nsec  pulses  and  is  estimated  to  be  ~  1.2  x  10-^  (W/cm^)-* 
at  steady  stated In  fact,  the  photothermal  effect  can  be  10  times 
larger  if  a  sample  has  10  times  larger  a.  For  laser  lines  that  fall  into 
the  absorption  band  of  a  mixture,  the  dn/dT  may  be  enhanced  due  to  a 
resonant  effect. 


III.  NONLINEAR  OPTICAL  MECHANISMS 


Nonlinear  self-focusing/defocusing  and  nonlinear  interface  will  be 
discussed  in  this  section.  Because  the  nonlinear  effect  of  an  AKM  is  about 
three  orders  of  magnitude  slower  than  the  photothermal  effect  at  the  same 
response  speed,  we  will  only  use  the  latter  as  an  example.  The  results  can 
be  qualitatively  applied  to  the  cases  involving  AKM,  as  long  as  the  time 
scale  is  increased  according  to  their  difference  in  response. 

A.  SELF-DEFOCUSING^8) 

When  a  thin  sample  of  liquid  mixture  is  exposed  to  a  pulsed  laser 
producing  a  TEMQO  laser  beam  with  a  radially  symmetric  Gaussian  intensity 
distribution  (Fig.  4),  the  sample  is  most  strongly  heated  at  the  center  of 
the  beam,  where  the  intensity  is  greatest,  and  consequently  forms  a  lens¬ 
like  optical  element.  For  most  materials,  the  increase  in  temperature 
lowers  the  refractive  index  so  that  the  optical  path  is  shorter  at  the  beam 
center,  i.e.,  it  is  equivalent  to  a  divergent  lens. 

As  shown  in  the  literature,  the  thermally  induced  effect  is  time 

dependent;  the  lens  requires  a  finite  time  to  reach  a  steady-state 

condition.  The  steady-state  focal  length,  f(»),  of  a  thermal  lens  produced 

2 

by  a  Gaussian  laser  beam  of  spot  size,  w^  ,  at  the  sample  has  been 
derived:^  ^ 

T  2 
irj  <w. 

f(oo)  P(dn/dT)a£  (30) 

where  J,  k,  P,  (dn/dT),  a,  and  l  are  defined  as  before.  This  expression 
assumes  that  all  of  the  abosrbed  radiation  is  converted  to  heat.  As  the 
thermal  lens  approaches  steady  state,  its  focal  length  f(t)  is 
characterized  ^2)  jjy  ^e  following  expression: 

f(t)  -  f(-)  (1  +  tc/2t)  (31) 
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Propagation,  2 


Figure  4.  Optical  configuration  for  calculating  the  focal  length  of  a 
thermal  lens  and  the  transmission  through  aperture  A. 
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where  tc  is  a  time  constant  given  by: 


t 

c 


“fop 

4< 


(32) 


and  p  and  Cp  are  defined  as  before.  Characteristic  time  constants  are  in 
the  range  of  tens  of  microseconds  to  several  seconds, ^ depending  on  the 
spot  size  of  the  laser  beam  and  the  thermal  properties  of  the  sample. 


Although  the  magnitude  of  tc  at  first  appears  hopelessly  long  for  an 
ns  switching  application,  the  time  needed  to  reduce  transmission  to  a  safe 
level  can  be  much  faster.  Using  the  optical  configuration  shown  in  Fig.  4, 
we  were  able  to  calculate  the  time-dependent  reduction  of  transmission 
through  an  aperture  A  for  a  1-W  laser  focusing  to  a  spot  20  pm  in  radius  on 
a  100- pm  sample.  The  sample  was  assumed  to  have  properties  similar  to 
benzene,  except  for  the  absorption  coefficient.  Also,  heat  transfer  from 
the  liquid  to  the  sample  cell  was  not  considered.  The  transmission  for  A  = 
0.002  and  0.055  cm  was  estimated  to  drop  to  <  50%  in  less  than  1  and  24  us, 
respectively  (Fig.  5). 


For  the  case  under  consideration,  i.e.,  t  «  tc,  f(t)  becomes 


f(t)  =  f(«)  tc/2t 


ttJ  w,  pc  1 

1  _B _ 

8P(dn/dT)  a£  t 


(33) 


i.e.,  the  response  is  proportional  to  w^  ,  p,  and  cp,  and  inversely 
proportional  to  P,  dn/dT,  a,  and  £.  Very  short  switching  times  are 
possible  by  properly  combining  these  parameters. 


B.  NONLINEAR  INTERFACE ^ 1 5) 


While  the  self-defocusing  mechanism  requires  a  focused  beam  to 
provide  lens-like  effect,  the  nonlinear-interface  effect  can  deal  with 
broad  beams.  Consider  an  infinite  plane  wave  incident  upon  a  plane 
boundary  separating  two  different  optical  media:  a  linear  medium  a  and  a 
nonlinear  medium  b  (Fig.  6).  In  addition  to  the  incident  wave,  there  will 
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Figure  5.  Dynamic  nonlinear  transmission  through  aperture  A  in  Fig.  4  due 
to  thermally  induced  self-defocusing  for  A  *  0.055  (upper 
curve)  and  0.002  (lower  curve)  cm.  Absorption  loss  in  the 
liquid  is  14* ;  overall  transoittances  in  the  absence  of  a  laser 
(i.e.,  t  -  0)  are  85  and  74*,  respectively. 


be  a  reflected  wave  and  a  transmitted  wave.  For  the  coordinates  shown  in 
Fig.  6,  the  interface  is  the  x-y  plane,  the  interface  normal  is  the  z  axis, 
and  the  wave  vectors  of  the  waves  lie  in  the  x-z  plane. 


The  reflected  angle  0r  and  the  transmitted  angle  ©t  are  related  to 
the  incident  angle  0^  as 


er  -  e4 


.  sin©. 

9t  *  Sln‘  (-T“) 


(34) 

(35) 


where 


a 


ng  and  n^  are  the  indices  of  refraction  of  the  two  media,  and  n  is  the 
relative  index  of  refraction. 

In  general,  EM  waves  can  be  resolved  into  two  components  with 
polarizations  perpendicular  to  each  other.  Reflectance  and  transmissivity 
at  the  interface  depend  on  the  polarization.  For  transverse  electric  (TE) 
polarization,  in  which  the  electric  vector  of  the  incident  wave  is  parallel 
to  the  interface,  the  transmissivity  and  the  reflectance  are  given  by 


TE 


1  -  R 


TE 


TE 


cos  0, 


(37) 


J  2 

-  /  n  - 


sin^  0. 


cos  0^  +  /  n^  -  sin“^0 


(38) 


*> 
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Transmitted  wave 


Figure  6. 


Schematic  diagram  of  the  nonlinear  interface.  Media  a  and  b 
occupy  the  regions  Z  >  0  abd  <  0,  respectively. 
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For  transverse  magnetic  (TM)  polarization,  in  which  the  magnetic  vector  of 
the  incident  wave  is  parallel  to  the  interface,  the  transmissivity  and  the 
reflectance  are  given  by 


a 


ft 


ft 


ft 


e 


c 


TM 


1  "  rtm 


(39) 


2  /  2  2 

n  cos©^  -  /  n  -  sin  9^ 


— 2 - 7  '!  '1  ~ 

n  cos 9^  +  /  n  -  sin  9^ 


(40) 


Transmissivity  and  reflectance  of  an  EM  wave  without  specific  polarization 
are  given  by  their  average  value 


T  ■  V2  *  Tm) 

(41) 

'■V2(RTE*RTM>- 

(42) 

Equations  (38)  and  (40)  clearly  show  that  transmissivity  and 
reflectance  depend  on  n  and  9j  for  both  the  TE  and  TM  cases.  Although  the 
actual  functions  are  different,  transmission  and  reflection  for  both 
polarizations  have  a  common  feature:  total  internal  reflection  occurs  for 
both  cases  when  n  <  1  and  9^  >  ©c,  where 

9  -  sin_1(n) .  (43) 

c 


Only  cases  in  which  medium  b  is  nonlinear  and  has  a  negative  n2« 
i.e.,  njj  ■  n^  +  n2l,  will  be  discussed  here.  However,  since  Eqs.  (38)  and 
(40)  depend  only  on  n,  the  results  car  be  easily  extended  to  cases  with 
positive  n2  as  well  as  to  cases  where  "a"  is  the  nonlinear  medium  and  "b" 
is  the  linear  medium. 

Substituting  n^,  into  Eq.  (43)  gives  the  power  density-dependent  9C 

as 
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(44) 
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9  ( I) 
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sin 


-1 


"bo  *  "21 
_  • 
n 

a 


In  other  words,  when  the  incident  power  density  increases,  the  critical 

angle  9C  decreases.  If  the  incident  angle  9^  is  selected  to  he  less  than 

9  (0)  ■  sin  *  (n,  /n  ) ,  a  low-power  beam  will  be  transmitted  through  the 
c  bo  a 

interface;  but  as  power  density  increases  and  the  critical  angle  decreases 
and  becomes  closer  to  9^,  transmissivity  is  reduced.  Further  increases  in 
I  shift  9C  to  angles  less  than  9j  and  result  in  total  internal  reflection, 
l.e.,  the  optical  density  of  the  limiter  becomes  infinite  since  no 
transmission  is  possible. 

Figure  7  shows  the  theoretical  transmission  through  a  nonlinear 
Interface  for  various  incident  angles.  At  large  9j,  transmissivity  changes 
are  more  abrupt  when  I  Increases.  For  example,  when  9^  *  80°  and  n^0/na  *  0.988, 
the  low  power  transmissivity  is  higher  than  80%.  But,  a  change  in  n  of  as 
little  as  0.002  due  to  an  intensity  increase  completely  blocks  the  high- 
power  Incident  beam.  At  9^  -  60°  and  n^x>/na  "  0.885,  transmissivity  can 
also  be  higher  than  80%.  However,  n  must  change  by  at  least  0.021  to 
produce  total  Internal  reflection  of  the  interface.  Therefore,  the  larger 
9^  is  preferred  for  optical  switching  because  the  threshold  power  density 
needed  to  cause  switching  from  high  transmissivity  to  total  internal 
reflection  is  smaller.  However,  a  selection  of  large  9^  will  limit  us  to  a 
small  field  of  view.  The  changes  in  indices  of  refraction  needed  at 
various  9j  to  switch  the  interface  from  80%  transmission  to  total 
reflection  for  a  TE  wave  are  listed  in  Table  1.  The  changes  required  for  a 
TM  wave  are  somewhat  smaller. 


Some  degree  of  switching  can  also  be  achieved  by  an  increase  in  0t, 
which  is  Induced  by  an  increase  in  the  intensity  dependent  nfl. 

Substituting  n^  into  Eq.  (35),  we  obtain 
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sin 


-1 


sin9t  \ 

TTno-J* 


(45) 
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Table  1 


Change  in  Indices  of  Refraction,  An,  Needed  for  Switching  an 
Interface  from  80%  Transmission  to  Total  Internal 
Reflection  for  a  TE  Wave 


0^  (degrees) 

n80%U) 

"c(b) 

An<c> 

89 

0.99987 

0.99985 

-0.0002 

• 

88 

0.99948 

0.99939 

-0.0009 

87 

0.99883 

0.99863 

-0.00020 

86 

0.99792 

0.99756 

-0.00036 

4t 

85 

0.99675 

0.99619 

-0.00056 

80 

0.98704 

0.98481 

-0.00213 

75 

0.97097 

0.96593 

-0.00904 

• 

70 

0.94873 

0.93969 

-0.00904 

65 

0.92057 

0.90631 

-0.01426 

• 

60 

0.88683 

0.86603 

-0.0208 

a)  Relative  index 

b)  Relative  index 

at  which  T  -  80% 

at  which  total  internal  reflection 

occurs 

c)  An 


nc  "  n80% 


e 


e 
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As  shown  in  Table  2  for  0^  *  60°,  the  resultant  0t  is  77.56°  for  a  low- 
power  beam.  A  change  in  n  of  only  0.0135  will  deflect  the  transmitted  beam 
5°,  so  that  becomes  82.56°.  With  0^  *  70°  and  0t  *  82.09°  for  a  low- 
power  beam,  the  required  change  in  n  can  be  even  smaller  (0.00782)  to 
deflect  the  beam  the  same  5°.  As  shown  in  Fig.  3,  such  a  change  in  n  can 
occur  in  100  usee  or  less  for  a  0.38-W  laser  focused  to  40  un  in 
diameter.  Therefore,  the  nonlinear-interface  phenomenon  is  also  fast. 


Table  2 

Change  in  Indices  of  Refraction,  Ajj,  Needed  for  the 
Transmitted  Beam  to  Bend  5° 


Low  Power 

High  Power 

Rj  (degrees) 

n80%  et 

0t-et+5°  ( A,.,) 

75 

0.97087 

84.16 

89.16 

-0.0049 

70 

0.94873 

82.09 

87.06 

-0.0078 

65 

0.92057 

79.90 

84.90 

-0.0107 

60 


0.88683 


77.56 


82.56 


-0.0135 


IV.  CONCLUSION 


We  have  analyzed  two  dynamic  nonlinear  optical  effects,  artificial 
Kerr  and  photothermal,  and  reviewed  two  nonlinear  phenomena,  self- 
defocusing  and  nonlinear  Interface.  The  generic  results  are  summarized 
below.  Please  see  Volume  II  for  their  device  applications. 

•  As  far  as  we  know,  we  are  the  first  to  derive  the  dynamic 
nonlinear  coefficient  n£  of  AKM. 

•  We  found  that  the  figure  of  merit  of  AKM  should  be  defined  as 
n2^ao^o  an<*  not  n2^ao  or^ginally  used  in  Ref.  6  because  aQ 
cannot  be  infinitely  increased. 

•  The  dynamic  nonlinear  coefficient  n^  is  larger  for  smaller 
spheres  under  the  constraint  of  the  same  scattering  loss. 

•  Although  both  AKM  and  absorptive  liquid  mixtures  are  highly 
nonlinear  optical  materials,  the  photothermal  effect  of  the 
latter  is  about  three  orders  of  magnitude  higher  under  the  same 
f ocus  condition. 

•  Both  the  self-defocusing  and  nonlinear-interface  phenomena 
produce  a  response  faster  than  100  usee  under  reasonable  power 
levels . 


t  r 
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